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The Differential Emission Measure of Dynamic Coronal Loops 
Spiro K. Antiochos 

Institute for Plasma Research, Stanford University 
Abstract 


Although most of the observations of soft x-ray and UV 
emission from stellar sources can be adequately explained in 
terms of static models of coronal loops, there are nov; several 
observations in soft x-rays of a solar-like flare from a stellar 
flare source, (e.g. Johnson 1981, Kahler 1981; Stern, et al. 

We discuss the effects of time dependant phenomena, such 
as flare energization and decay, on the temperature and density 
structure of the transition region and, in particular, on the 
form of the differentia] em.ission measure. We find that unlike 
the case of the static models, the form of the differential 
emission measure can be used to determine the important physical 
mechanisms in the dynamic models. 

I. Introduction 


Since the advent of solar x-ray observations with high 
spatial resolution, (in particular the data from Skylab) , there 
has been considerable interest in the so-called "quasi -static” 
model of coronal loops, (e.g. Landini and Fossi 1975, Rosner et. 
al., 1978, Craig et. al. , 1978, Hood and Priest 1979, Vesecky et 
al., 1979). The important assumptions in the "quasi-static" 
model are that (a) the plasma in static and (b) the magnetic 
field dominates the plasma so that a one-dimensional treatment 
(i.e. parallel to the field) is sufficient. Additional assump- 
tions that are usually made are that; (c) the conductive flux 
is small compared to the saturated value so that the collisional 
thermal conductivity is valid (Spitzer 1962) ; (d) all radiation 

is optically then so that the losses calculated by e.g. Cox and 
Tucker (1969) or Raymond and Smith (1977) , can be used; and (e) 
the coronal heating can be expressed as a simple function of 
density and temperature, (usually two power laws) . 

The key result of the static model is that the plasma tem- 
perature and density profiles, along the loop, are essentially 
independant of the assumed form for the heating function. The 
quantity obtained from the observations is the differential emis 
Sion measure, which we define as: 

- 1 



where A is the cross-sectional area of the loop, n is the 

electron density, and ^ In T is the temperature scale height 

d s 

along the loop. The static model predicts that: 

T ® with 6 » 3/4 + jl/2 , (2) 

where 8. measures the slope of the radiative loss coefficient 
(see Cox and Tucker 1969) , i.e. 

A(T) ~ t'^ (3) 


Since ^ ~ 1/2 for transition region temperature ( Rosner et . 
al., 1978), ve find that 6=1, for a static loop. 

• 

Note that, as stated above, the differential emission mea- 
sure is independent of the heating function parameters, eqtn. (2). 
Hence, observations of C are ineffective for determining the 
nature of the coronal heating mechanism; however, they do pro- 
vide a stringent test of the assumption that the coronal plasma 
. is static. If a value for 6 is observed that is not compatible 
with eqtn. (2) , then this argues strongly in favour for a dynamic 
model. In practice, eqtn, (2) provides only an upper limit for 
6 since the observed emission is unlikely to be due tc» a single 
loop, and the effect of adding the contributions from a collec- 
tion of loops is to decrease the value cf 6 , (e.g. .Antiochos 

1980) . 

The observed value of 5 for solar quiet and active regions 
are generally i 1 , and hence are compatible with the static 
model, (Noyes, et. al., 1970, Raymond and Doyle 1981). Unfor- 
tunately, there is little data in the EUV region on stellar 
coronae, so that 6 has not been accurately determined for stars 
other than our sun. What data is available also appears to be 
compatible with eqtn. (2) (Zolcinski, et . al . , 1981). Solar 
flares, on the other hand, are generally not compatible with 
eqtn. (2). Sevecal authors (Dere, et. al. , 1977, Underwood et. 
al., 1978, Widing and Spicer 1980) have reported a large value 

for 6, - 3, in the temperature range, 10^ - T ~ 1C'. This 
value cannot be reconciled with the static model even if a col- 
lection of loops is observed. 

Of course, flares are not static so that it is not surpris- 
ing that their values for 6 do not agree with static models. 

The important point here is that observation of the form of the 


differential emission measure can be used to infer the presence 
of dynamic effects, even if the observations do not have suffi- 
cient temporal or spatial resolution to observe these effects 
directly. 

In the next section we derive the value for 6 " that is 
predicted by the dynamic models. It turns out that this value 
is sensitive to the particular physical mechanism that dominates 
the dynamics. 

II . Dynamic Models 

The dynamics that we consider are assumed to be due basic- 
ally to a time dependant coronal heating function. This time 
dependence could represent either an impulsive event such as a 
flare, or a more gradual evolution such as is expected to occur 
during the formation and decay of an active region loop. Note 
that we will not consider extremely rapid events such as chromo- 
spheric evaporation by a beam of energetic electrons, (the so- 
called ''primary” evaporation) . For such processes, the form of 
the differential emission measure will clearly depend on the 
details of the energization, (e.g. the exact time profile, energy 
spectrum, etc. of the electron beam) , and will be different for 
every particular event. In addition, non-thermal effects, such 
as departures from ionization equilibrium, are likely to be 
important so that the simple differential emission measure as 
defined in (1) is no longer relevant to the observations. For- 
tunately such processes, (e.g. the hard x-ray bursts in solar 
flares) , are typically short-lived compared to the total flare 
duration as measured, for example, by the soft x-ray bursts in 
solar flares. 


We identify three important phases for the "thermal" evolu- 
tion of a flare-like loop, that predict a characLeristic signa- 
ture for the form of the differential emission measure, i.e. for 
6. These phases are characterized by the ratio, R (T) , of the 
radiative flux from loop plasma at temperature T, to the con- 
ductive flux at that temperature: 


R 


( T) = 


n^ A (T) 
lO"^ 


T 


(4) 


where 


is the temperature scale height, as given in eqtn.(l). 


We define the three phases as: (a) conduction dominated, for 

which R << 1 throughout the loop, i.e. for T 10^; (la) con- 
densation dominated, for which R << 1 in the corona, i.e. for 
T > 10°, but R ^ 1 in the transition region, T < 10®; and (c) 
radiation dominated, for which R - 1 throughout the loop. 



IIA Conduction Dominated Models 


In the early stages of a flare event, we expect conduction 
to dominate the dynamics. This results from the following argu- 
ment. Assume a coronal loop is initially static. In this case 
R • 1 throughout the loop (Vesecky, Antiochos and llnderwood 
1979) , so that radiation and conduction are approximately equal 
throughout. If the energy input to the loop begins to rise 
rapidly, as in a flare event, then the initial response of the 
plasma is to increase its temperature. The density cannot rise 
until material is evaporated up from the chromosphere. Since _ 
the conductive flux increases rapidly with temperature, P « T' , 
whereas the radiative losses do not, (in fact, they decrease 
somewhat for 10^ $ T ^ 10^) , the effect of the increased energy 
input is to enhance conduction over radiation, i.e. R decreases. 
If the increase in energy input is significant, as in a flare, 
then R becomes << 1 throv jut the loop and conduction domin- 
ates. 


This results in a large heat flux entering the chromosphere 
from above. Since radiation is insufficient to dirsipate this 
flux, it must result in a large enthalpy flux, and hence .mass 
flux, being driven up the loop. This is the process of "chromo- 
spheric evaporation” as described by Antiochos and Sturrock 
(1978) . Under some simplifying assumptions, the temperature 
gradients, and therefore the differential emission measure, can 
be calculated quite readily for the evaporative model. We 
assume a one-dimensional model, neglect the effect of gravity, 
and assume that the heat flux is well below its saturated value 
so that the motions are highly subsonic. In this case the 
relevant equations are: 




J_P = 0 , ( 6 ) 

8s 


3/2 4P 


1 .5 

A 8s '2 


ApV - lO'^AT^/^ ^'^) = e(t) (7) 

C & 


where p is the mass density, A(s) is the cross-sectional area 
of the loop and s(t) is the energy input function, assumed to be 
a function of time only. Note that we have not included the 
radiative losses in the heat equation since these are presumed 


to be small compared to the other terras in the equation. 

Imposing the boundary conditions that at the loop foot- 
points the velocity and heat flux become small, eqtn. (7) 
implies that: 


2 ~ ds 


( 0 ) 


Substituting (8) into (5) yields a nonlinear equation for T(t,s) 
that can be solved by separation of variables. Investigating 
the form of the solutions of this equation, it turns out that 
over most of the temperature range, (i.e. except for tempera- 
tures very near the maximum temperature in the loop) , the heat 
flux, 


F 




dS 


(9) 


Since 



(10) 


1/2 

the evaporative model predicts that f, « T over most of the 
temperature range. Note that in deriving (10) v.’e have used the 
fact that the loop is isobaric, eqtn. (6), and that from (9), 


F « 


( 11 ) 


In addition, we have used the result that the size scale for 
variations of the magnetic field, and A(s) , are large compared 
to the size scale of the transition region. 

Our main conclusion is that conduction driven mass motions 
tend to decrease the slop of C » (5 = 1/2) , from the value given 
by the static models (6 * 1). Unfortunately, the effect of ob- 
serving a collection of static loops is also to decrease the 
value of 6 . Hence, the observation that 6 = 1/2 does not 
necessarily imply that conduction driven evaporation is occurr- 
ing, unless the observations also indicate that the loop density 
is increasing with time. 


and 


pv = const 


• i 


( 12 ) 

(13) 


pv 


2 


+ 


P 


const . , 


_d_ 

ds 


(2 pv 




„ 5/2 

ds^ 


n^A(T) , 


(14) 


whore we have assumed that the loop cross-section varies insig- 
nificantly over the transition region. 

The condensation model described by equations (12) -(14) has 
been investigated by Antiochos and SturrocJc (1982). The impor- 
tant results of their analysis are that (a) the coronal plasma, 

7 

T ~ 10 K, cools primarily by mass motion, rather than by conduc- 
tion or radiation; and (b) in the region where the flow is super- 
sonic, 6 = 3.5, The latter result is due primarily to the fact 
that the flow is isochoric rather than isobaric. 


The )cey aspect of the condensational model, for stellar 
purposes, is that it predicts a larger value for 6 than can be 
explained by the static models. In addition, we expect that it 
is valid during most of the decay phase of a flare, because the 
time scale, t, for cooling by mass motions increases with de- 
creasing T. In particular, » 

T oc i “ ^“1/2 
C 

where C is the sound speed of the coronal plasma. Hence, if 
stellar flare emission is detectei, at all, it is li)cely to be 
that predicted by the condensation dominated model. 

III. Radiation Dominated Models 


In the later stages of flare decay, we expect that radia- 
tion cooling dominates throughout the loop. This can be seen by 
comparing the ratio, R' , of the radiative cooling rate to that 
by mass motion: 


R' 


n'A (T) 

pc 


n/T' 


(16) 


Since, in the condensation phase, most of the coronal energy 
loss is due to mass motion, rather than conduction or radia- 
tion (Antiochos and Sturrock 1982), the coronal plasma cools 
primarily by adiabatic expansion so that 


n 


5/3. 


hence, R’ “ T 


- 1/2 


(17). 


This result implies that eventually radiation will dominate even 



in the corona. However, since R* has a weak dependence on T, 
radiation will not dominate until the coronal plasma has cooled 
to quite low temperatures, probably < 10 ®, and has lost the 
bulk of its energy. Hence, the total amount of emission during 
the radiation dominated phase is likely to be negligible com- 
pared tc that radiated during the condensation phase. 

In any case, we have calculated the differential emission 
measure to be expected for a radiative dominated loop (Antiochos 
1980) , and it turns out to predict a value for 6 of only ~ 1 
This value is not sufficiently different from that given by the 
static models so that it can be used as an observational test. 

Of more importance ■ for the radiation dominated model is the 
rapid decay expected for the emission. 

IV. Conclusions 


In summary, we find that the form of the differential emis- 
sion measure is a sensitive indicator of the important physical 
processes in the dynamic models. We tabulate the slope of C (T) 
predicted by the various models below: 


Model 6 


Static 

1.0 

Conduction 

0.5 

Condensation 

3.5 

Radiation 

1.5 


The model that has most significance for stellar observations is 
the condensation dominated model which predicts a large value 
for 5 and also large supersonic velocities. 
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IIB Condensation Dominated Models 


Once the flare energy input ceases, the temperature in the 
loop will begin to drop, whereas the density still increases as 
a result of conduction driven evaporation. These two effects 
act to increase the ratio R(T) . Eventually this ratio reaches 
unity at some point in the loop, and chromospheric evaporation 
ceases since now radiation is sufficient to dissipate the down- 
ward heat flux. The temperature at which R first reaches unity 

is ^ lO^K as can be seen from Fig. 1, where we indicate the be- 
havior of R{T) for a typical conduction dominated loop. 



Figure 1. R(T) for the conduction dominated models. 

When R(T) becomes of order unity at T - 10^ , the loop 
plasma at this temperature begins to cool very rapidly because 
the downard heat flux is no longer sufficient to sustain the 
large radiative losses. As a result of this rapid cooling, a 

pressure difference between the material at 10^ and that at -lO^ 
is created which, in turn, generates large downward velocities. 
This process may be thought of as the opposite of chromospheric 
evaporation in that the loop plasma cools and ondenses onto the 
chromosphere . 

Assuming, again, a one-dimensional model and, in addition, 
assuming that in the transition region where the velocites are 

large (supersonic), 10^ 5 T 5 10^*^, the flow is approximately 
in a steady-state, the relevant equations become; 


